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' The first fermion family might play a special role in understanding the physics 

of flavour. This possibility is suggested by the observation that the up-down splitting 
within quark families increases with the family number: m u ~ m^, m c > m s , Wit ^> 
m&. We construct a model that realizes this feature of the spectrum in a natural way. 
The inter-family hierarchy is first generated by radiative phenomena in a sector of heavy 
isosinglet fermions and then transferred to quarks by means of a universal seesaw. A 
crucial role is played by left-right parity and up-down isotopic symmetry. No family 
symmetry is introduced. The model implies m^/m^ > 0.5 and the Cabibbo angle is 

forced to be ~ \J md/m 8 . The top quark is naturally in the 100 GeV range, but not 
too heavy: m t < 150 GeV. 

Inspired by the mass matrices obtained in the model for quarks, we suggest an 
ansatz also including charged leptons. The differences between u-, d- and e-type 
fermions are simply parametrized by three complex coefficients e u , €d and e e . Ad- 
ditional consistent predictions are obtained: m s =100-150 MeV and m u /md < 0.75. 
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1. Introduction 



One of the major issues in modern particle physics is to understand the origin of the ob- 
served pattern of fermion masses and mixing. In the Standard Model all the flavour struc- 
ture is determined by uncalculable Yukawa couplings. A hypothetical Theory of Flavour 
should allow a calculation of these couplings, or somehow constrain them. When thinking 
of such a theory, one should keep in mind the following qualitative features: 

i) An inter-family mass hierarchy which is stronger for the up quarks than for the down 
quarks [1]. The plot in Fig. 1 suggests for the i-th family quark masses 

~ el^mt mf ~ e^mf (1) 

where e" 1 = 200 - 300 and e^ 1 = 20 - 30. 

ii) The CKM mixing matrix V is close to the unit matrix, and it has a hierarchy in its 
off-diagonal entries: V u b ~ e d , V c b ~ £d and V us ~ e 1 / 2 . 

The inter-family hierarchy makes the radiative picture of fermion mass generation 
[2-4] attractive. In particular, the quark mass pattern (1) can be interpreted in terms of a 
charge diagonal radiative cascade. In this view e u and are the loop expansion parameters 
respectively in the up and the down sector [4] . Leptons, on the other hand, seem to evade 
the simple rule rrii ~ e 3 ~ l ms : and this suggests that their masses may be more difficult to 
understand (see Fig. 1). We therefore focus first on the quark sector. 

When discussing quarks, we can also exploit experimental information on the CKM 
matrix. The relation V ~ 1 implies that the Standard Model Yukawa matrices T u and 
T d are somehow aligned. This property might hint that some "isotopic" symmetry in- 
terchanging u- and (i-quarks has played a role in the mass generation phenomenon. (A 
horizontal family symmetry could also be responsible for the alignment of T's. However 
the idea of an up-down symmetry is simpler.) Let us suppose that such a symmetry indeed 
exists, and let us call it I u a- In the low energy theory, i.e. the Standard Model, I u d plays 
the role of the usual custodial symmetry [5]. In the limit of exact I u a it would be Y u = Fj, 
so that we would have V = 1 and the up and down sectors unsplit. Within the radiative 
picture the breaking of I u d implies the appearance of small mixing angles V u b ~ e d and 
Vch ~ Q. 1 ) 

As far as the spectrum is concerned, Fig. 1 shows that I u a is more badly broken for 
the heavier families. By focusing on the first family, we don't see a big violation of I uc i: 



1 ^ ) However, also V us has to be 0(ed)- This is a generic problem of the direct radiative 
scenarios [3-4]. As it was shown in Ref. [6], the correct value V us ~ 0.22 implies a choice 
of parameters which spoils the validity of perturbation theory. 
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m-a/md = 0(1). What we consider interesting is the case in which I ud becomes exact at 
very high energies, namely the case of spontaneous or soft breakdown. With this picture 
in mind, it is suggestive to think that the masses of u and d are somehow related to an 
energy scale Ai at which I u d is still good, while the masses of the second and third families 
are respectively related to lower scales A 2 , A 3 at which I ud is no longer as good. Thus, the 
first family seems to play a special role in mass generation. This is directly opposite to the 
common radiative picture [2-4], in which the third family (i.e. the heaviest) is the starting 
point. However, by taking the top and bottom masses as the seeds of mass generation, it 
is difficult to understand their large splitting. 2 -* On the other hand, if we could start from 
u and d, we would have better chance to understand t-b splitting. Suppose that u and d 
are indeed the starting point, and that an expression like eq. (1) holds for l/m^ instead of 
rrii, namely 

y<' d -^/m. (2) 

Then we have m u ~ ma ~ m and nit/nib ~ (e^/e^) 2 ~ (ni c /m s ) 2 ^> 1. In this way, 
the splitting between up and down quark masses in Fig. 1 is understood by means of one 
parameter (e<f/e u ) > 1. We call the above formula for l/m^ the inverse hierarchy pattern. 

In this paper we explore the idea of inverse hierarchy for fermion masses. In Section 
2, we discuss an illustrative model for the quark sector. The model is a particular case of 
the class of models discussed in Ref. [7]. These models naturally avoid the key problems 
of the previous models of radiative mass generation (see footnotes 1 ' 2 -'). A set of isosinglet 
heavy quarks, Q-fermions, in a one to one correspondence with the ordinary ones (<?'s), 
is introduced. The mass matrices of the q's M u and Md are induced by a seesaw mixing 
with the Q's [8]. In terms of the mass matrices M UjD of the Q's we have essentially the 
inverse proportionality 3 -' M u oc M^ 1 , Md oc M^ 1 . Thus the smallest masses m u and md 
are indeed related to the highest scale: m" 1 oc M v and m^ 1 oc M D , where M v and M D 
are the largest eigenvalues respectively in M v and M D . These matrices are generated by 
a radiative mechanism analogous to that in Ref. [4]: the largest eigenvalues M UjD appear 
at tree level, while the other ones are induced by charge diagonal loop effects. Moreover, 

2 ) In Ref. [3-4] the large t — b splitting is introduced at tree level by hand. A difficulty 
of charge diagonal radiative models is also the appearance of flavour-changing neutral 
currents (FCNC), generally contradicting the experimental bounds [4]. 

3 > The seesaw mechanism [9], universally extended to quarks and charged leptons was 
also explored in the papers in Ref. [10,11]. The inverse hierarchy, however, corresponds to 
the spirit of the original paper [8], where it was first suggested, though in the context of 
horizontal symmetry. 
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the relation M v ~ M D arises naturally as a consequence of an up-down symmetry. Then 
m u ~ rridi so that the pattern of eq. (2) is realized. 

In the model we discuss below, the field content and the I u a symmetry imply strictly 
Mu = M D . The model also has Plr and CP symmetries. These, together with I u d, 
are softly (or spontaneously) broken in the scalar potential. The effects of this breaking 
are reflected in the fermion sector via loop corrections only. As a consequence, the mass 
matrices are more constrained than those in Ref. [7]. In particular, the splitting between u 
and d is related to the enhancement of the Cabibbo angle above its natural radiative value 
0(ed)- The consequences, including formulae for the top quark mass and the Cabibbo 
angle, are discussed in Section 3. 

The quark mass matrices obtained in the model are a very attractive realization of 
the inverse hierarchy. Inspired by the form of those matrices, in Section 4 we simply 
postulate an inverse hierarchy ansatz, which includes charged leptons. It turns out that 
the first family (u, d, e) really plays the role of a mass unification point, with its splittings 
understood by the same mechanism that enhances V us . This ansatz provides remarkable 
predictions for the light quark masses, consistent with the present knowledge on these 
quantities. 

Finally, in Section 5 we discuss our results. 



2. A model 

Let us consider the simple left-right symmetric model, based on the gauge group 
G LR = [SU(2) L x U(1) L ] x [SU(2) R x U(1) R ] x £/(%_ r , suggested in [7]. The left- 
and right-handed components of usual quarks qi = (ui, di) and their heavy partners Qi = 
Ui, Di are taken in the following representations 

(1/2,0,0,0,1/3) q Ri ~ (0,0, 1/2, 0,1/3) 

(0,1,0,0,1/3) U Ri -(0,0,0,1, 1/3) (3) 

(0, -1, 0, 0, 1/3) D Ri ~ (0, 0, 0, -1, 1/3) 

where i = 1,2,3 is the family index (the indices of the colour SU(3) C are omitted). We 
also introduce an additional family of quarks 4 ) 

p L ~ (0, -1/2, 0, 3/2, 1/3) p R ~ (0, 3/2, 0, -1/2, 1/3) 

( 4 ) 

n L ~ (0, 1/2, 0, -3/2, 1/3) n R ~ (0, -3/2, 0, 1/2, 1/3). 

4 ) Notice that, by adding the obvious lepton fields (with B — L = —1) to the quarks of 
eqs. (3-4), the theory is free of gauge anomalies. 



q L i ~ 
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The scalar sector of the theory is given by 





(1/2,0,0,1,0) 


H R ~ 


(0,1,1/2,0,0) 




(0,-2,0,0,-2/3) 


T U R ~ 


(0,0,0,-2,-2/3) 


TdL ~ 


(0,2,0,0,-2/3) 




(0,0,0,2,-2/3) 


(j) ~ 


(0,1/2,0,-1/2,0) 




(0,2,0,-2,0) 


o u ~ 


(0,1/2,0,1/2,-1) 




(0,-1/2,0,-1/2,-1) 



(5) 



where the T-scalars are colour triplets or sextets. We impose CP, Plr and I u a discrete 
symmetries. The left-right symmetry Plr [12], which is essentially parity, and CP act in 
the usual way. The "up-down" symmetry I ud is defined by 



U l ,r D LjR , Pl,r n LjR , R — > A1 R , 

h l , r ^h l , r , t? jR »tZ jR , 0^0*, n u ^n d 



(6) 



where A^ )R are the gauge bosons of U(1)l,r- Then the most general Yukawa interactions, 
consistent with gauge invariance, I u d, Plr and CP are 

d = T ij {q Li U Rj H L + q Ll D RJ H L + (L <-> P)} + h.c. 

£ 2 = ^-{T UL U Li CU Lj + T dL D Ll CD LJ + (L <-> P)} + h.c. (7) 
£3 = {hi((f)*U Li p R + <f)D Ll n R ) + h($*p L p R + m L n R ) + (L <-> P)} + h.c. 

where C is the charge conjugation matrix. The coupling constants /i, ftj, A^, (z,j = 
1,2,3) are real, due to CP-invariance. In what follows we do not make any particular 
assumption on their flavour structure. We only suppose that they are typically 0(1), just 
like the gauge coupling constants. Notice that A is antisymmetric (symmetric) when the 
T-scalars are colour triplets (sextets). Without loss of generality, by a suitable redefinition 
of the fermion basis, we can always take 

h 2 ,3 = 0, Ai 3 = 0, Fij = if % < j . (8) 

We assume that the discrete symmetries CP, Plr and I ud are softly broken by bilinear 
and trilinear terms in the Higgs potential 5 ) . The trilinear terms are given by 

U = KT: l T ur § + A d T2 L T dR Q* + h.c. (9) 



5 ) Actually, these symmetries can be broken spontaneously by suitable scalars [7]. 
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where A Utd are mass dimensional coupling constants. We suppose that these couplings are 
complex, and also that A u ^ A d . In this way, CP, Plr and I ud are broken. 

The v.e.v.'s ($) = vq>, ((f)) = and (Q u ) break the U(1) L x U(1) R x U(1)^_ T 
symmetry down to the usual U(1)b-l- B — L = Y l + Y r + B — L. (Notice that we must 
have (Qd) = 0.) Then, interactions £3 in eq. (7) give rise to mass terms for p and n, 
and also for the Q-fermions of the "first family" U\ and D\. It is interesting to consider 
the limiting case v$> ^> v^. In this case, the tree- level mass eigenvalues are given by 
the seesaw formula: M p — M n — hv^ and M Ul ^ Dl ecause of the 

relation M <C M p n , the mass matrices obtained after including radiative effects depend, 
in general, on a smaller number of parameters. In fact, we shall see below that a soft 
hierarchy hv^/hiv^ ~ 3 is already enough for our purposes. 

The coloured scalars T UL ,T UR and T d L,T d R get mixed due to interaction terms in eq. 

(9) . At this point the interactions in £ 2 trigger the radiative mass generation chains 
U\ — > U2 — > Us and Di — > D 2 — > -D3 (see Fig. 2). The two-loop corrected mass matrices 
for the Q-fermions can be written as 

M UjD =M{ A + e u4 X T P 1 X + C u4 \e u4 \ 2 A t2 A A 2 

. (10) 

+ tu, d \ T (^, d A T2 + < d r T r) A a + e ., d A T A {^ d \ 2 + < d r T r) a} 

where A = diag(l, 0, 0) is the tree- level term. In terms of the mixing angle $, Ujd and mass 
eigenvalues M^' d , M"' d of the scalars 7^' d -T^' d , the expansion parameters are 

^ 16yr (11) 
C u , d =-(1 + H(r u , d )) . 

Here cu U) d are the arguments of the mixing terms A Ujd in eq. (9), while H is a real function 
given in the Appendix. Ct is a colour factor equal to 1/2 for a triplet T and to 1 for a 
sextet T. The first three addenda in eq. (10) represent respectively the contributions of 
graphs (a), (b) and (c) of Fig. 2. They are rank 1 matrices determining the eigenvalues 
of M to be O(M), 0(eM) and 0(e 2 M). The terms proportional to (3 Ujd and a U}d in eq. 

(10) arise from the graphs in Figs. 2d-e. Apart from pieces which are formally 3-loop, i.e. 
(9(ect 2 , e/3 2 , ect/3), these terms simply determine a multiplicative redefinition of the 1-loop 
contribution (Fig. 2b). (We remind that ct's and /3's are loop coefficients ~ 1/16/t 2 .) 
This means that in 2-loop order their effect on the eigenvalues is only multiplicative, i.e. 
Mi — > Mj(l + 0(a, 0)) [3]. For our purposes, these terms are therefore negligible. 
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Some comments on eqs. (10-11) are in order. Eq. (11) is valid in the regime M 2 <C 

2 2 

M"' , M"' <C M 2 . In this regime the loop-diagrams of Fig. 2b-c are dominated by loop 
momenta between M™' d and M"' d . Thus, the result is completely determined by the ratios 
r UjC i, and the predictivity is increased. The choice of the intermediate range for M± d is 
also motivated by the fact that both for M^ d < M and for M^ d > M p the loop integrals 
are suppressed. In the first case the suppression factor is given by ~ (M+'7M) 2 , in the 
second by ~ (M p /M"' ) 2 . (In fact M p acts cut-off of the mass insertions in Fig. 2 
[4].) Notice also that (M/M p ) 2 ~ (hiv^/hv^) 4 <C 1 is already implied by our assumption 
> v^. For instance, hiv$/hv$, ~ 3 already gives (M/M p ) 2 ~ 10 -2 . 

The function H(r) = H(l/r) has the maximum at r = 1 and goes to zero mono- 
tonically when r — > oo (see Appendix). However, in a reasonable range (1 < r < 10) it 
is H(r) ~ 0.3 with good accuracy. Thus, in what follows we shall set C u ^ = C ~ 0.65. 
Notice that the two loop term proportional to C u ^ in eq. (10) is real. This is because the 
mass insertions on the two scalar lines in Fig. 2c are complex conjugate of each other. 

As a final comment, notice that, in the basis of eq. (8), we have 

/O(l) 0(6) 0(e 2 )\ 
M oc O(e) O(e) 0(e 2 ) (12) 
\0(e 2 ) 0(e 2 ) 0(e 2 )/ 

Then, in lowest order, the eigenvalues are given by the diagonal entries of M, and their 
ratios are 0(1) : 0(e) : 0(e 2 ). 



3. The quark mass matrices 

We now discuss the masses of the ordinary quarks q. The v.e.v's (H L ) = (0, v L ) and 
(H L ) = (0,v R ), with v R » v L = 174 GeV, break the intermediate SU(2) L x SU(2) R x 
U(1)b-l symmetry down to U(l) em , and at the same time mix the g's to the Q's via the 
interactions Ci in eq. (7). The mass matrix of the -u-type quarks is written in block form 

as 

and analogously for the <i-type ones. When M UyD ^> v R , the resulting mass matrix M U) d 
for the ordinary quarks is given by the seesaw formula [8-10] 

M u , d = v L v R TM-} D T T (14a) 
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In this way the inverse proportionality of eq. (2) is realized [7,8]. The seesaw limit M UjD 3> 
v R is certainly very good for light quarks, since their masses must be much smaller than 
v L . However, since m t = 0(v L ), we expect the mass of its Q-partner M T to be of the order 
of v R (remember that the Y^ are considered to be 0(1)). Thus, to evaluate nit, we need 
the mass matrix without the restriction M v ^> v R . This is given by 

{MujM^r^v^Y-^v^M.JiYY^M,^ + l}^)" 1 . (146) 

Notice that, when v R 3> M v ,m this equation gives the obvious result M Uj d = v L Y. On the 
other hand, when v R <C M v ,o, eq. (14b) reduces to the seesaw formula. 

The v.e.v. v R ~ M T plays the role of the Flavour scale Ap. Below the scale v R the 
effective theory is the minimal standard model. In fact v R TM~* D Y T = M Ujd /v L are the 
standard model Yukawa matrices Y Ujd arising at the decoupling of the heavy Q states. 
The Q-fermion masses are such that M T <C M B <C ... <C M v = M D = M. Then the ratio 
M/v R is given by v L /m u ~ 10 5 . Obviously, FCNC's are well controlled in such a theory 
[8] due to natural flavour conservation within the standard model [13]. Flavour changing 
effects are suppressed by inverse powers of the Q masses. Phenomenologically, the most 
stringent constraints come from K° — K° mixing and K L — > n + fi~ . It turns out that the 
typical scale of suppression for these processes is given by M. This is readily seen in basis 
(8): the light quarks u,d mix only with the heaviest Q-fermions U, D. Consistency with 
the data is obtained for M > 100 TeV. Since M/v R ~ 10 5 , we have that v R ~ 1 - 10 TeV 
is already consistent with the experimental data. Therefore, it is in principle possible to 
have the W R bosons and the lightest Q-fermion, T, in the TeV range, with implications at 
SSC/LHC. 

It is useful to first study M u ^ d in the seesaw limit (14a). The exact formula (14b) will 
only be relevant to evaluate mt. Once again, the inverse matrices are easier to analyse. 
Neglecting a and (3 terms, as explained above, eq. (10) gives 

(M^)" 1 = m" 1 ! A + e u , d 7 T A 7 + C\e u4 \ 2 7 t2 A 7 2 + • • } (15) 

where 7 = TAr -1 and m = T^VlVr/M. We also used the relation (r T ) _1 Ar _1 = 
A, true in the basis of eq. (8). Notice that, in general, 7 is no longer symmetric or 
antisymmetric. However Tr7 = Tr A, Det7=DetA and, because of eq. (8), 713 = A13 = 0. 
Eq. (15) is then written in the same form as eq. (12), which is easily diagonalized to lowest 
order in e. 

There is, however, a subtlety involving the Cabibbo angle. Let us consider eq. (15) 
restricted to the first two families at 0(e). We have 
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This gives, in lowest order, sin6>c — Cd7ii7i2 and d/s ~ £d1i2 2 - These relations imply 
e d 7n 2 ~ sin 2 $c(s/d) ~ 1, which looks like trouble for perturbation theory [6-7]. How- 
ever, 7n = Ai2r 2 i/r 2 2 is not a coupling constant: its value can be enhanced by the ratio 
r2i/r22, without trouble for perturbativity [7]. This is an advantage of the seesaw mech- 
anism: the "sandwiching" between T's in eq. (14a) allows the mass matrix in eq. (15) to 
deviate a little bit from the rigid radiative form of eq. (12). In particular, when r 2 i is 
larger than Fn and 1?22, the 1-2 entry gets enhanced. This effect is reflected in eq. (15-16) 
by the fact that 711 is large. Notice that, by supposing e^A^ 2 ~ d/s = 1/20, we only need 

r 2 i/r 22 = 4-5. 

To accommodate a large Cabibbo angle, the matrix (15) must be diagonalized con- 
sidering that ed7n 2 = 0(1). With an obvious notation, we indicate by -u, <i, ... the quark 
mass eigenvalues at the flavour scale Ap. Then we get 



(m/u) = |1 + e u 7n : 
im/d) = |1 + tdlii 



(m/c) = 
(m/s) = 



IM712' 



1 + e«7n 2 | 



M7i~ 2 



1 + ed7n^ 



(m/t) = CK| 2 7 i2 2 723 2 
(m/b) = C|e d | 2 7i2 2 723 2 



and the Cabibbo angle is 



V us = 711712 



m' 



7ii V \ e d\ 



$u,d = arg(l + e* d 7n 2 )- 



^ e i{u u +8 u -uj d -8 d ) 
C 



(17b) 



By dropping terms proportional to e u ^n 2 in the above equations, we would get the naive 
lowest order result. (Notice that we have kept also e u 7n 2 terms. However, we shall see that 
£u «C €d, so that those terms are not relevant). Notice also that the third family masses 
are simply given by the 3-3 entry in eq. (15), i.e. they do not get relevant contributions 
from off-diagonal terms. This is because their mixings V c b and V u b must be of the natural 
radiative size, namely V c b ~ ~ d/s, s/b and V u b ~ e 2 , ~ d/b. 

The following comment is in order. In contrast to the quantities in eqs. (17), V c b and 
V u b get lowest order contributions which are not expressed in terms of e. One of these 
comes from Fig. 2d-e: these graphs contribute in lowest order to V c b and V u b, but not 
to the quantities in eq. (17). In particular, V c b gets contributions 0(a,(3), whereas V u t, 
gets 0(ect, e(3). This is consistent with the radiative picture of mixings, since a, (3 are loop 
parameters ~ e. 

Another contribution comes from the renormalization of T's. These effects are given 
at 1-loop by triangle graphs involving three Yukawa vertices in C\ in eq. (7). The mass 
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splittings between up and down Q-fermions determine this renormalization to differ by a 
finite amount in the two sectors. This is an additional source of mixing. It is remarkable, 
however, that a contribution appears at 1-loop only in the mixing between the second and 
third family, namely in V c b. This is a consequence of the fact that the first family (u, d) is 
only connected to the first heavy family (U, D), which is unsplit in lowest order. Then the 
vertex graphs involving u and d are equal and do not contribute to the mixing in lowest 
order: the relation M v = M D protects V us and V u \, from 1-loop corrections and thereby 
maintains V u b two loop effect. 

The above observation selects the present model from the wide class discussed in Ref. 
[7], not only as a predictive one (see below), but also as the only one in which the mixing 
angles V c b and V u b follow the natural radiative pattern: V c b ~ 1-loop, V u b ~ 2-loop. 6 ) 

Let us now discuss the implications of eqs. (17). We have seen that V us is enhanced by 
seesaw effects above the radiative value 0{td)- Eqs. (17) show the connection between the 
large value of the Cabibbo angle and the u-d splitting. Neglecting terms of order \e u /ed\, 
we get 



!d 

i/ ~ 

S 



1 - -e li 
d 



(18) 



l li 

Thus, without any hypothesis on 5d, we have V us in the right range: \V US \ ~ e d . For 
instance, for u/d < 0.6 and s/d = 20 [1], we obtain 0.14 < \V US \ < 0.26. Notice that the 
experimental value of |V^ S | requires large values 8d ~ 1. On the other hand, Sd contributes 
appreciably to the CP phase in the CKM matrix, even though there are other independent 
contributions. These are related to the effects on V u b and V c b that we discussed above. 
However, barring accidental cancellations, we can conclude that the model favours a large 
CP-phase, in accord with experiments [14]. 

Let us now discuss the mass eigenvalues. From the masses of the first two families in 
eq. (17a) we get 

1^1 = —. (19) 

e d cu 

On the other hand b/t = \e u /ed\ 2 , so that we have a mass formula 

HI) 2 - < 20 > 



6 ) The relation M v = M D does not generally hold in the models of Ref. [7] . It was con- 
sidered to be broken at tree level. Consequently, in order to have a small enough V u b, some 
Yukawa coupling in C\ (namely I^i) must be taken about one order of magnitude smaller 
then the others. In Ref. [7] the role of the renormalization of T's had been overlooked. 
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This equation, valid in the seesaw limit of eq. (15), is an appealing expression for the 
top mass. However, even when exact, eq. (20) does not translate into a sharp top mass 
prediction, because of our poor experimental knowledge of R = (uc/ds) 2 . By using the 
correct mass matrix of eq. (14), eq. (20) is reduced to an inequality 

(i/A t ) 2 = {i/x b R) 2 + (i/r 33 ) 2 > {i/x b R) 2 (21) 

where A^ = mt^/v L are the standard model Yukawa couplings. Eq. (21) is valid at the 
"Flavour scale" v R , and, to discuss its implications, the running of masses needs to be 
considered. For our purposes, we only need to consider the combined effects of QCD and 
of the top Yukawa coupling [15] (electroweak gauge couplings modify our conclusions by 
less than 5%). In what follows we assume a s (M z ) = 0.11 and use the estimate of the QCD 
running of masses given in ref. [16] at 2-loop accuracy. There are two implications from 
eq. (21). 

i) An upper bound on the top mass . This value is determined by the maximal allowed 
values for R and r 33 . Taking m u /md < 0.8 and m c /m s < 10.7 from the values given in Ref. 
[1,17] we have R < -R max = 74 (Notice that R is not affected by the running of masses, 
when electroweak couplings are neglected [15].) On the other hand, T 33 is bounded by 
requiring the perturbativity of the model: T 33 < 0.5 — 1. This upper bound requires 
some explanation. From the definition of 7, and from the basis choice eq. (8), we have 7 -* 
T2i/r 33 = (ed7n72 3 )/(edAi2A2 3 ). Then, from the expressions for d/s and sin6>c in eq. 
(17), we have €^711 ~ 1, while from the b-mass we have e^72 3 2 ~ 0.06. On the other 
hand, e^Ai2A2 3 is a true loop expansion parameter in eq. (10), so that we demand that it be 
smaller than 0.1 — 0.05 ~ d/s. Putting all these constraints together we get r2i/r 33 > 2 — 4. 
Finally, by demanding the reasonable perturbative upper bound T21 < 2 we obtain the 
bound r 33 <, 0.5 - 1.0. 

The results are shown in Table la for different values of A F . Notice that for = 
10 12 GeV the heaviest states of the model lay just below the Planck mass. However, for 
the sake of demonstration we also kept A F = 10 16 GeV. As we see the top quark cannot 
reasonably be heavier than 150 GeV. For instance, for A F = 10 8 GeV, in order to push m t 
to 170 GeV we should take T 33 ~ 2.2. Such a value, as discussed above, would lead to a 
large T 2 i > 4 and the loss of perturbative control. In fact, at the same value of A F: even 
in the limiting seesaw case, i.e. T 33 — > 00 in eq. (21), the top quark cannot weigh more 
than 176 GeV. 

7 ) Actually this relation is true only in the case of a triplet T. However, the consequences 
in the case of a sextet do not differ appreciably. 
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ii) A lower bound on \^R . This is a consequence of the CDF bound m t > 91 GeV [18] 
(see Table lb). Notice that the lower bound on \fR = (m u m c /mdm s ) is given, in practice, 
by R > [mt / rribimt)} ^ 30 = -R m i n . Then, by using the maximal value c/s = 10.7 [1], we 
get 

u/d > VR min (s/c) > 0.5 (22) 

The determination of u/d in chiral perturbation theory is still controversial, essentially 
because the axial U(1) A is broken by the anomaly and instanton effects may mimic a 
substantial portion of the u quark mass [19]. The range of u/d allowed in our model is 
consistent with the values reported by Leutwyler in Ref. [17] . On the other hand, Donoghue 
and Wyler [20] give the estimate u/d = 0.3 ± 0.1. If this last estimate turns out to be the 
correct one, then our model is ruled out, unless c/s J> 14. By taking m c (l GeV) < 1.4 GeV 
[1], this bound translates into m s <, 100 MeV, which seems too small a value for m s . (Even 
though it is very difficult to determine the absolute scale of the strange quark mass [1]). 

4. Inverse Hierarchy Ansatz with charged leptons 

Lepton masses have a mixed behaviour, m M /m e ~ e u and m T /m M ~ e^, which seems 
hard to explain by means of only one expansion parameter (see Fig. 1). On the other hand, 
we saw that the large value of the Cabibbo angle requires a deviation from the genuine 
radiative pattern in the quark sector. In fact the off-diagonal entries in the mass matrix 
affect sizably the masses of the first two families. Here we want to show that the same 
mechanism allows to explain charged lepton masses by means of one complex expansion 
parameter e e . 

We simply postulate the inverse hierarchy form of eq. (15) for all mass matrices 8 ) 

{M k )~ l = m- 1 {P 1 + e k 7 T A 7 + M 2 7 t2 A 7 2 } (23) 

where k = u,d, e. This ansatz is simple: the difference between -u-, d- and e-type fermions 
is parametrized by the e k only. The first family plays the role of a mass unification point. 

The eigenvalues of quarks at the Flavour scale are the same as those in eq. (17) 
(with C = 1) while for leptons we have the analogous formulae 

(m/e) = |1 + e e7 n 2 | (m//i) = J^ 712 ' (m/r) = \e e \ 2 ll2 2 723 2 . (24) 

I J- + e e 7n I 

8 ) We set C = 1 by rescaling 723. This does not affect the generality of the ansatz. 
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Comparing the above equation with eq. (17a), we see that the experimental relation r ~ b 
implies |e<f| — \e e \. Then |e e |7n 2 ~ 1 as in the down sector, so that e and fx receive in 
general sizable corrections from this term. In particular, when the phase of e e is close to 
zero, we have that fi is enhanced and e is decreased as required by the experimental data 
(Fig. 1). Eq. (24) gives a mass formula involving leptons and analogous to eq. (20) for 
quarks. We can write it in terms of down-type quark and lepton masses 




Running down this relation from the Flavour scale Ap, we get a prediction for m s md 

m s m d = m e m M /^y\V /2 ^ /2 - 277 z^rfy '* ( MeV) 2 . (26) 

Where r] s is the QCD running of masses from m t = 130 GeV to 1 GeV, and r\h is the 
analogous quantity from m t to m^. The remaining z t and r\ F account for the running from 
A F to m t : 7] F = (a s (m t ) / a s (A F )) 4 / 7 represents the effects of pure QCD. While z t gives 
the effects of the top quark Yukawa coupling [15]. For m t < 180 GeV the effects of z t in 
eq. (26) are always < 5%, so that we neglect it. The quantity m s md for different values 
of hp and the corresponding results for the light quark masses are displayed in Table 2. 
We have let A F go up to 10 16 GeV. 9 ) The values reported correspond to a s (M z ) = 0.11, 
and the r/'s have been deduced from Ref. [16]. Notice that the prediction for m s (and 
hence ma) is increased by about 30% when a s (M z ) = 0.13. The lower bounds on m u /md 
correspond to the limit R = R m in = 30. We also get an upper bound for m u /md- Indeed, 
from eqs. (17a) and (24) we obtain 

li.^lH^.Wl^) 1 ' 2 © 1 ^^^!^)^) 172 (27, 

where x = {r\hr\Fm T j mb) 1 ^ = 0.9 — 1.1, for Ap varying between 10 4 GeV and 10 16 
GeV. To derive the inequality in eq. (27) we have used eq. (25) for \e e /ed\- Then rather 
independently on we obtain u/d = |1 + ed7n 2 |/|l + e u 7n 2 | < 0.75, by considering that 
|ed|7n 2 ~ \V US \ 2 (s / d) ~ 1 due to eq. (17b) and \e u /ed\ < 1/6 due to eq. (19). 

The ansatz (23) gives correct predictions by taking the first family as a mass unification 
point (namely, by taking the same normalizing factor 1/m for u, d, and e type fermions). 

9 - ) Notice that even in a model like the one of the previous section it is meaningful to 
run up to such a scale. The fact is that, in comparing quarks to leptons, QCD running 
becomes effective above v R . 
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We can also reverse the argument, by testing different factors 1/M, in eq. (23), instead 
of the same 1/m. It is then interesting to see what values for Mj's are implied by our 
experimental knowledge on fermion masses. We have (Md/M u ) 3 = (t/b)(ds/uc) 2 and 
(Md/M e ) 3 = (r/b)(ds/e/j,) 2 . For instance for Ap = 10 12 , after rescaling masses, we get 

Md -i/3 ( mt 4.25GeV \ 1/3 / m s 1.35 GeV \ 2/3 f 06 ™d\ 2/3 
M u ~ Zt Vl20GeV m b J VlSOMeV m c J \'m u ) 



Ma -i/6 / 4.25GeV \ 1/3 / m s \4/3 /W 2/3 
M e ~ Zt \ m b J U30MeV/ \ m 

where we have put the lepton masses to their values. It is remarkable to see that, by 
varying the masses of eq. (28) in their allowed experimental ranges, all Mj's remain equal 
within a factor of 2. 



5. Discussion 

Motivated by the observation of the inverse hierarchy pattern l/m™' d ~ e lTd/ m ^ n 
the quark spectrum and by the fact that an up-down symmetry could be responsible for 
the smallness of the CKM mixing angles, we have constructed a realistic model for quark 
masses. The key features of this model are: 

i) A universal seesaw mechanism [8,10]: the masses of the ordinary quarks q are induced 
via their mixing with ultraheavy families of weak isosinglet fermions Q. 

ii) The mass matrices of the Q's are generated by a charge diagonal radiative cascade 
[4,7]: the I s * family gets its mass at tree level, while the 2 nd and 3 rd get their masses 
at 1- and 2-loops, respectively. 

iii) The left-right P LR and up-down I u d discrete symmetries. These are softly or 
spontaneously broken in the Higgs potential. 

The heaviest Q-family (U, D) is unsplit because of the I u d symmetry. The lighter Q's 
get their masses from loop-diagrams involving the breaking of I u d and are thereby split. We 
considered a particular, but reasonable range for the masses of the scalar fields involved in 
mass generation (see the discussion below eq. (11)). In this range, the Q mass matrices are 
simply determined by two expansion parameters e u and e^. These are in general complex 

10 - ) In fact, we used a discrete up-down symmetry just for simplicity. I u d can be eas- 
ily extended to a local "isotopic" invariance via the embedding U(1)l x U(1)r x I u d C 
SU(2) < 1 x SU{2) ( ^ t . Again, for simplicity we have imposed CP-invariance. Our results are 
unchanged in the case of complex Yukawa coupling constants. 
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and e u ^ ej, due to the breaking of CP, P LR and I uc i- The SU(2) R breaking scale plays 
the role of the Flavour scale. The Q masses are such that v R ~ M Us < M D3 < ... < M. 
The ratio v R /M is determined by m u /v L ~ 10~ 5 . As a result, the q : s mass matrices, given 
by a seesaw mixing with the Q's, have the inverse hierarchy form displayed in eq. (15). 

Once again we would like to stress that in our approach the ordinary light quarks 
are just spectators of the phenomena that determine the flavour structure. This structure 
arises in a sector of heavy fermions and is transferred to the light ones at their decoupling. 
This explains why flavour changing effects are suppressed in low energy phenomena. 

The mass matrices of eq. (15) reproduce the quark spectrum (see Fig. 1) and mixing 
angles in a very economical way. The big differences between u- and <i-type quark masses 
are essentially determined by one parameter |e^/e u | = ~ 6 — 8. Upper bounds for 
rrit and md/m u are then obtained. Rather independently on the Flavour scale = 
v R , by invoking reasonable perturbativity requirements, we have rrit < 150 GeV. On 
the other hand, without further assumptions, we obtain (m u m c /mdm s ) > 6. Then the 
relation m u /md > 0.5 is obtained by using the strange quark mass range given in Ref. 
[1]. The experimental input m u /rrid < 1 forces the the Cabibbo angle to be large: V us ~ 
(md/rris) 1 / 2 instead of V us ~ md/m s as one generally expects within the radiative scenario. 
Moreover, the correct value V us ~ 0.22 does not imply the loss of perturbativity, in contrast 
to the previous models of radiative mass generation. 

The model shows that the inverse hierarchy pattern can be obtained in an ordinary 
field theory. Having in mind the idea of some symmetry between quarks and leptons, 11 ^ eq. 
(15) can be generalized to charged leptons, in the inverse hierarchy ansatz of eq. (23). In 
this ansatz the differences between u-, d- and e-type fermions are simply parametrized by 
three complex numbers e u , and e e . In fact, in eq. (23) we have kept the same notation of 
eq. (15) just for convention: actually A, A = 7 t At and A = 7 t2 At 2 parametrize the 
most general set of three symmetric and real matrices of rank 1. Then, eq. (23) contains 
11 relevant parameters to calculate 13 observables (9 masses + 4 angles). Hence, we have 
the two relations in eq. (20,25). (Notice, however, that even by taking A hermitean these 
relations remain; this simply means that CP-invariance is not essential for our results.) 
The predictions can be phrased as m s = 100 — 150 MeV and m u /md = 0.35 — 0.75. 

n ) The question of whether the S77(4) quark-lepton symmetry of Pati-Salam [21] could 
provide this ansatz requires further speculations. The T-scalars of SU(4) contain lepto- 
quarks, that violate in general the charge diagonality of the radiative cascade. Also, 577(4) 
or any GUT extension brings neutrinos into the game. The related physical consequences 
are under study. 
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We find it amusing that the radiative idea, implemented with the simplest symmetries 
(P, CP, "isotopic" and "lepton-quark" symmetries), can explain the key features of the 
fermion mass spectrum and weak mixing. Notice, that we did not exploit any horizontal 
structure, in contrast to all known predictive frameworks for fermion masses (see e.g. 
[11,22]). Clearly, a clever horizontal structure would only enhance the predictive power of 
our approach. 

Last but not least, we wish to remark that our approach can automatically solve the 
strong CP-problem, without introducing an axion. In fact, due to P/CP-symmetries the 
coefficient Oqcd of GG is zero. In addition, the contribution to O, that arises from the 
determinant of the whole mass matrix of the fermions q, Q and p, n, is vanishing at tree 
level. This is a consequence of P/CP and of the seesaw structure [23]. Higher order 
corrections to the mass matrix can provide a 6-term in the range 10~ 9 — 10~ 10 which is 
of interest for the search of the the neutron electric dipole moment. 
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Appendix 

We report the expression of the function H{r). Let us define: 

f°° dx 1 

f(r)= — -ln(x + l) = -{lnrln II- r| +3?Li 2 (r)} Al 

Jq x(x + r) r 

where 3?Li2 is the real part of the Euler dilogarithm. Then H is given by: 

H(r) = -^[f(r) + f(l/r)-^-]. A.2 
in r o 

H is a positive function which reaches its maximum at r = 1: H(l) = n 2 /3 — 3 ~ 0.290. 
In the range 1 < r < oo, H decreases monotonically and very slowly to zero (see Fig. 3). 
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Figure captions 

Fig. 1 Logarithmic plot of fermion masses as a function of the family number. Points corre- 
sponding to fermions with the same electric charge are joined. Leptons are given by 
the dashed line. The value m t = 130 GeV has been assumed. The masses correspond 
to a renormalization scale \x = 10 6 GeV. 

Fig. 2 a) tree level diagram giving mass to U via mixing with p. b) one-loop correction to 
the mass matrix of the Q-fermions. c,d,e) two-loop corrections. The analogues of d) 
and e), where a T R and a H L are respectively exchanged, are not displayed. 

Fig. 3 Plot of H{r). 
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